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^ : Abstract 
O 

Csl ■ We derive upper limits on the ratio /grb/ccsnI^) = -Rgrb(^;)/^ccsn(2;) = /grb/ccsn(0)(1 + 2;)"*, 



the ratio of the rate, i?GRB) of long-duration Gamma Ray Bursts (GRBs) to the rate, -RccsN, of 
core-collapse supernovae (CCSNe) in the Universe {z being the cosmological redshift and a > 0), 
by using the upper limit on the diffuse TeV-PeV neutrino background given by the AMANDA-II 



1- ^ • experiment in the South Pole, under the assumption that GRBs are sources of TeV-PeV neutrinos 

Q , produced from decay of charged pions produced in interaction of protons accelerated to ultrahigh 

H ■ 

tj^ ' energies at internal shocks within GRB jets. For the assumed "concordance model" of cosmic star 

formation rate, RsF, with i?ccsN(-z) oc Rsf{z), our conservative upper limits are /grb/ccsn(0) < 
^ . 5.0 X 10"'^ for a = 0, and /grb/ccsn(0) < 1.1 x 10^^ for a = 2, for example. These limits are 

I already comparable to (and, for a > 1, already more restrictive than) the current upper limit on this 

10 ' ratio inferred from other astronomical considerations, thus providing a useful independent probe 

of and constraint on the CCSN-GRB connection. Non-detection of a diffuse TeV-PeV neutrino 
, background by the up-coming IceCube detector in the South pole after three years of operation, 

^ for example, will bring down the upper limit on /grb/ccsn(0) to below few xlO^^ level, while a 

^ • detection will confirm the hypothesis of proton acceleration to ultrahigh energies in GRBs and will 

potentially also yield the true rate of occurrence of these events in the Universe. 
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I. INTRODUCTION 



Detection of supernova (SN) features in the afterglow spectra of several long duration 
(typically > 2 s) Gamma Ray Bursts (GRBs) in the past one decade has provided strong 
support to the hypothesis that a significant fraction, if not all, of the long duration GRBs 
arise from collapse of massive stars; see, e.g., Refs. [2,y] for recent reviews. The observed 
SN features in the GRB afterglow spectra are similar to those usually associated with core- 
collapse supernovae (CCSNe) of Type Ib/c (see, e.g., 0, Q]). The total energy (corrected 
for beaming) in keV-MeV gamma rays emitted by typical long-duration GRBs is of order 
Iq5i gj,g^ which is roughly the same as the total explosion energy seen in typical CCSNe, 
although there exists considerable diversity in the energetics of both the SN and the GRB 
components in the SN-GRB associations observed so far. In particular, the estimated ex- 
plosion energies of the SNe associated with the GRBs observed so far seem to be somewhat 
larger than those of normal SNe, leading to this "special" class of SNe being sometimes 
referred to as "hypernovae" . 

The broad class of observational results on SN-GRB associations can be understood 
within the context of the "coUapsar" model in terms of a simple phenomenological pic- 
ture (see, e.g., jol) in which the core-collapse of a massive Wolf-Rayet star gives rise to 
two kinds of outflows emanating from the central regions inside the collapsed star: (a) a 
narrowly collimated and highly relativistic jet that is responsible for the GRB activity, the 
jet being driven, for example, by a rapidly rotating and accreting black hole formed at the 
center in the core-collapse process, and (b) a more wide-angled, quasi-spherical and non- 
relativistic (or at best sub-relativistic) outflow that goes to blow up the star and gives rise 
to the supernova. The energies channeled into these two components may in general vary 
independently, which may explain the diversity of energetics in the observed SN-GRB as- 
sociations. Actually, depending on the energy contained in it the "GRB jet" may or may 
not be able to penetrate through the stellar material and emerge outside. Indeed, the fact 
that the SN-GRB associations observed so far involve CCSNe of Type Ib/c, but not of Type 
II, may be due to the inability of the GRB-causing jet to penetrate through the relatively 
larger amount of outer stellar material in the case of Type II SN as compared to that in SNe 
of Type Ib/c [7]. Considering various factors that may govern the energy channeled into the 
GRB-causing jet, such as the mass and rotation rate of the black hole, accretion efficiency. 
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efficiency of conversion of accretion energy into collimated relativistic outflow, and so on, 
Woosley and Zhang gI have obtained a rough lower limit of ~ 10^^ erg/ s for the power 
required for the jet to be able to emerge from the star. This is consistent with the energetics 
of the GRB components of the SN-GRB associations observed so far. 

While SN-GRB associations strongly support the stellar core-collapse origin of most long- 
duration GRBs, clearly, not all core-collapse events may result in a GRB — the latter 
depends on whether or not the core-collapse event actually results in a "central engine" 
(a rotating black hole fed by an accretion disk in the above mentioned phenomenological 
picture, for example) that is capable of driving the required collimated relativistic outflow. 
In other words, while every long-duration GRB would be expected to be accompanied by a 
core-collapse supernova the reverse is not true in general. 

What fraction of all stellar core-collapse events in the universe produce GRBs? Methods 
based on astronomical observations generally indicate the ratio between the cosmic GRB 
rate and the cosmic Type Ib/c SN rate, /cRB/SNibc, to be in the range ~ 10^^ - 10^^ for 
a wide variety of different assumptions on various relevant parameters such as those that 
characterize the cosmic star formation rate (SFR), initial mass function (IMF) of stars, 
masses of Type Ib/c SN progenitors, the luminosity function of GRBs, the beaming factor 
of GRBs (associated with the fact that individual GRB emissions are highly non-isotropic 
and confined to narrowly collimated jets covering only a small fraction of the sky), and 
so on; see, for example, [ol, [l^ and references therein. The dominant uncertainty in the 
estimate of /cRB/SNibc comes from the uncertainties in the estimates of the local GRB rate 
and the average GRB beaming factor. However, irrespective of the exact value of the ratio 
/cRB/SNibc, it is clear that this ratio is significantly less than unity. This indicates that, 
apart from just being sufficiently massive stars, the GRB progenitors may need to satisfy 
additional special conditions. For example, it has been suggested [2] that the degree of 
rotation of the central iron core of the collapsing star and the metalicity of the progenitor 
star may play crucial roles in producing a GRB. 

In this paper, we discuss an alternative probe of the cosmic GRB rate that uses the pre- 
dicted high energy (TeV-PeV) diffuse neutrino background produced by GRBs and the ex- 
perimental upper limit on high energy diffuse neutrino background given by the AMANDA- 



II experiment in the South Pole 



Existence of a high (TeV-PeV) energy diffuse GRB 



neutrino background (DGRBNuB) due to pj interactions of (ultra)high energy protons ac- 
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celerated within GRB sources is a generic prediction 12] in most currently popular models 



of GRBs. This DGRBNuB is subject to being probed by the currently operating and up 



coming large volume (kilometer scale) neutrino detectors such as IceCube [13| , ANITA IJ] , 
ANTARES [15], for example. Since neutrinos, unlike electromagnetic radiation, can travel 
un-hindered from the furthest cosmological distances, the DGRBNuB automatically includes 
the contributions from all GRBs in the Universe. Thus, an analysis of the DGRBNuB is 
likely to provide a good picture of the true rate of occurrence of these events in the Universe. 
Indeed, as we show in this paper, the upper limits on /grb/ccsn(0), the ratio of the local 
(i.e., redshift z = 0) GRB to CCSN rates, derived here from the consideration of DGRB- 
NuB, are, for a wide range of values of the relevant parameters, already more restrictive 



than the current \m 



per limit on this ratio (~ 2.5 x 10 ^) inferred from other astronomical 



considerations [9|, [lOl]. Further, non-detection of a diffuse TeV-PeV neutrino background by 
the up-coming IceCube detector [isl in the South Pole after three years of operation, for 
example, will imply upper limits on /grb/ccsn(0) at the level of few xlO"^, while a detec- 
tion of the DGRBNuB will provide strong support to the hypothesis of proton acceleration 
to ultrahigh energies within GRB jets. 

Our use of the DGRBNuB in constraining the cosmic GRB rate is in the same spirit as 
efforts to constrain the cosmic star formation rate (and thereby the cosmic CCSN rate) by 
using the experimental upper limit (set by the Super-Kamiokande (SK) detector) |l^ on the 



predicted [17| low 
for example, Refs. 



few MeV) energy Diffuse Supernova Neutrino Background (DSNuB); see. 



18 
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201 1 . Now that the cosmic SFR including its absolute normalization 



and thereby the cosmic CCSNe rate have got reasonably well determined by the recent high 
quality data from a variety of astronomical observations (see, e.g., {2^) (which, by the way, 
predicts a DSNuB flux that is close to the SK upper limit, implying that the DSNuB is 
probably close to being detected in the near future), one can begin to think of using this 
SFR to constrain the ratio of the cosmic GRB rate to CCSNe rate by using the predicted 
DGRBNuB flux together with the recent upper limits on the diffuse high energy neutrino 
flux from neutrino telescopes. 

We should emphasize here that the upper limits derived in this paper actually refer to 
the ratio of the rate of GRBs to that of all CCSNe including those of Type Ib/c and Type 
II, although SN-GRB associations observed so far involve SNe of Type Ib/c only. It is 
known, however, that Type II SNe probably constitute as much as ~ 75% of all CCSNe; 
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see, e.g., 2l|. Thus, one can get the constraint on the GRB-to-SNIb/c ratio from the GRB- 
to-CCSNe ratio we obtain here by multiplying the latter by a factor of ~ 4. Conversely, for 
later comparison, we shall take the "observed" value of the ratio /grb/ccsn(0) to be in the 
range 2.5 x (lO"^ - IQ-^) [qI, 1^. 

Below, we first briefly review the calculation of the DGRBNuB spectrum in section [TTl 
The resulting upper limits on /grb/ccsn obtained by comparing the DGRBNuB with the 
current upper limit from AMANDA-II experiment are discussed in section IIIII for various 
values of some of the relevant GRB parameters. Finally, in section IIVI we summarize the 
main results and conclude. 

II. DIFFUSE HIGH ENERGY NEUTRINOS FROM GAMMA RAY BURSTS 



Starting with the original calculations of Waxman and Bahcall [12], the production of 
TeV-PeV neutrinos is widely accepted as a generic prediction of the fireball model of GRBs, 
provided, of course, that protons (in addition to electrons) are accelerated to ultrahigh 
energies within GRB jets. Reviews of the basic method of calculation of the expected 



neutrino flux from GRBs can be found, e.g., in 



22 



24 



25|. 



23| . Recent calculations of the GRB 



neutrino spectra can be found, for example, in 

For a given cosmological rate of occurrence of GRBs, the DGRBNuB flux can be calcu- 
lated by simply convoluting the neutrino production spectrum of individual GRBs with the 
GRB rate density as a function of redshift, integrating over redshift up to some maximum 
redshift, and averaging over the intrinsic GRB parameters. In this paper we closely follow 
the recent calculation of the DGRBNuB spectrum described in Ref. 25| with appropriate 
modifications for a possible enhanced evolution of the cosmic GRB rate in redshift relative 
to the cosmic SFR as indicated by a recent analysis of the Swift GRB data [2^. 

In the standard jet fireball model of GRBs (see, e.g., Refs. [22, 1271] for reviews), the fun- 
damental source of the observed radiation from GRB is the dissipation of kinetic energy of 
ultra-relativistic (Lorentz factor F ~ few 100) bulk flow of matter (caused by ejection from 
a "central engine") through formation of shocks which accelerate particles (electrons and 
protons) to ultra-relativistic energies. The shocks can form either inside the flow material 
itself due to collision of different shells of matter moving with different Lorentz factors ( "in- 
ternal shocks" ) or due to collision of the flow material with an external medium ( "external 
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shocks"). The emission of the observed prompt 7-rays from a GRB source is attributed pri- 
marily to synchrotron radiation (with possible additional contribution from Inverse Compton 
scattering) of high energy electrons accelerated in the internal shocks. 

It is expected that along with electrons, protons would also be accelerated at the internal 
shocks. Since synchrotron energy loss of protons is a slow process, protons can be accelerated 



to much higher energies than electrons. Indeed, it has been suggested 28|| that protons 



may be accelerated to ultra-high energies in GRB internal shocks and that these UHE 



protons may explain the observed ultra-high energy (UHE) cosmic rays (UHECR) 29|] with 
energies up to ~ 10^° eV. These UHE protons interacting with the photons within the GRB 
jet would produce high energy charged pions through the photo-pion production process, 
p + 7 n + 7i~^, and the subsequent decay of each charged pion would give rise to three 
high energy neutrinos (a u^, a and a z/g): 7r+ ^ /i+ + z/^ , yU+ — > e"*" + z/^ + z/g. 

The dominant contribution to photo-pion production comes from the A resonance, p + 
7 ^ A+ ^ n + 7r+, at which the jry interaction cross section peaks with a value 
10"^*^ cm^. This A resonance occurs at the proton threshold energy e^^j^ (as measured in the 
GRB wind rest frame — the "comoving frame" hereafter), which satisfies epth^-y ~ 0-3 GeV^, 
where is the comoving frame energy of the colliding photon. In the rest frame of the GRB 
source (i.e., the central engine), the above threshold condition is ep^th^^ ~ 0.3 GeV^, where 
r is the bulk Lorentz factor of the GRB wind. In each p7 interaction the pion takes away 
on average a fraction ~ 20% of the energy of the proton, so each neutrino from the decay of 
the pion carries ~ 5% of the energy of the initial proton, assuming that the four final state 
leptons share the energy of the decaying pion equally. Thus, for a typical photon energy 

~ 1 MeV and F = 300, say, we have e^^th ~ 3 x 10^ GeV, which will give rise to neutrinos 
of energy ~ 1.5 PeV. 

The observed prompt 7 ray spectra of most GRBs are consistent with photon spectra 



which are well described by a broken power-law |30|, ImI]: dn^/de^ oc e~^, with [3 ~ 1.0 for 
< e^b, and {3 ~ 2.25 for > e^h- For typical GRBs, the break energy e^b ~ 1 MeV. The 
normalized photon spectrum in the source rest frame (SRF) can be written as 

de-f It .1-25^ -2.25 



£76 £7 for £7 ^ £76 



where is the total photon energy density in the SRF. Note that quantities in the SRF are 
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related to those in the comoving frame (denoted by primes) by the appropriate powers of the 
Lorentz T factor. Thus, for example, e-y = Te'^, and f/y = f/^ = L^/Aixrl c, where is the 
photon luminosity in SRF, and r^^ = T'^ct^ is the characteristic "dissipation" radius where 
internal shocks are formed and from where most of the radiation is emitted, ~ (10~^ - 
10^^ sec) being the typical variability timescale of the emitted radiation. Note further that 
the quantities observed at earth (denoted by the superscript or subscript 'ob') are related 
to those in the SRF through appropriate powers of the redshift factor (1 + z). Thus, for 
example, = e,^/(l + z). 

We shall assume that at the internal shock protons are accelerated to a differential spec- 
trum, dup/dep cx Cp^. The total internal energy in the system, ftotai, is assumed to be 
distributed among electrons, protons and magnetic field as £e = ^e^totah = ^p^totai and 
£b = 'C-B^totai, respectively, with + C,p + C,b = 1- We further assume that electrons are effi- 
cient radiators, so that £e ^ = ^e^totai, where S-y = L^T^ is the total isotropic-equivalent 
energy of the emitted gamma ray photons, being the total duration of the burst. 

It is worthwhile noting here that in the fireball model the kinetic energy of the initial 
bulk fiow of matter is predominantly carried by protons, they being ~ 2000 times more 
massive than electrons. This kinetic energy then is converted into internal energy at the 
shock, whereby the energy is now shared by protons, electrons and magnetic field. The 
mechanism by which the energy, which is initially carried mainly in the form of protons, 
gets transferred to electrons (and magnetic field) is not clear, but the phenomenology of the 
observed radiation from GRBs requires a significant fraction of the total internal energy to 



be eventually carried by electrons (see, e.g., 23]). If this energy transfer from protons to 
electrons is very efficient, it may lead to equipartition of energy between them, i.e., C,p = ^e, 
but in general one may expect that C,p/C,e ^ 1- 

Now, with the proton and photon spectra specified as above, the photo-pion production 
interactions of the protons with the photons given by the spectrum in eq. ([T]) can be shown 



to give rise to the neutrino spectrum 
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.«.^x0.56x/.(e4.J; ^ (2) 

^^•^ ^ (e^/e^*)-2 for > e^. 



where fn{^p), the fractional energy loss of a proton to pions during the dynamical expansion 
time scale of the wind [l^, is to be evaluated at ep = 206^. For the photon spectrum given 
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by equation ([T]), fn{^p) has the form [25| 



0.88(ep/ep;,)i-25 for < e^^ , 
A(ep) = /o S (3) 
1 for ep > , 



with /o = 0.09L^,5i/(r^oot^_3e^6,Mev)- Here L^,5i = L^/(10^^ erg s"^), t^_3 = ti,/(10"^ s), 
r3oo = r/300, and e^^^MeV = e^^/ MeV. 

In equation ([2]) the factor ~ .^pf^totai = £p is just the internal energy contained in 
protons, of which a fraction goes to pions. The factor | comes from the fact that in 
P7 interactions vr+'s and vr^'s are produced with roughly equal probability and the three 
neutrinos from the decay of 7r+ together carry 3/4-th of the pion's energy. Finally, the 
factor 0.56 is an overall normalization factor. 

The spectrum ([2]) has two breaks: The first break at e^^ = 0.05 e^f, is caused by the break 
in /^(ep) at epb with 

ep, = 1.3 X lO^r^oo (e^b.Mcv)"' GeV, (4) 
which, in turn, is due to the break in the photon spectrum ([T]) at e^b- 



The second break is at e,^^, with [251] 



e,, = 2.56 X 10'ej'^B'^'L-l,{' T^oo tv,-3 GeV , (5) 

which is due to muon cooling; for neutrino energy above e^* the corresponding energy of the 
parent muon (coming from the decay of the pion) would be high enough that the character- 
istic timescale of its energy loss through synchrotron radiation ("cooling") would be shorter 



than its decay time scale. Following [25|] we shall assume = ^s, in which case e^^, becomes 
independent of these two parameters. For a given (which is an observationally measur- 
able quantity), the neutrino spectrum ([2]) then depends on ^p and ^e, but only through their 
ratio, ^p/^e, which we shall take to be a free parameter in our calculations below. 

With the neutrino spectrum from individual GRBs (in the GRB source rest frame) given 
by equation (j2]), the diffuse neutrino flux from all GRBs in the Universe, DGRBNuB, can 
be calculated as follows: 

Let drii.^e'^) denote the present number density of neutrinos with energy between and 
e"^ + (ie°^, which were emitted with energies between e^, and + dey from GRBs at redshifts 
between z and z + dz. Denoting by Rgkb{z) the GRB rate per comoving volume at redshift 
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2;, we can write 



dn^ie:^) = Rokb{z){1 + zf (^^Jz^ de. (1 + z^' . (6) 

In this equation the factor (l+z)^ on the right hand side converts the GRB rate per comoving 
volume to the rate per physical volume while the factor [1 + z)~^ accounts for the dilution 
of the number density of the produced neutrinos due to expansion of the Universe. 
Using the standard Friedmann relation 

- = - [Hoii + z)^n^{i + zY + nj,\ (7) 

(we shall use the standard ACDM cosmology parameters, Qm = 0.3 , Q\ = 0.7 and Hq = 
70 km/ s/ Mpc), the total differential flux of neutrinos, $(e°^), giving the number of neutrinos 
(of all flavors) crossing per unit area per unit time per unit energy per unit solid angle, due 
to all GRBs in the Universe up to a maximum redshift z^ax can be written as 

^ " ' ~ 47r def 

where '^^^j^''^ is given by equation with = {1 + z)e^. We assume that, because of 
the long cosmological baseline, neutrino flavor oscillation distributes the original neutrinos 
equally into all three flavors. 

The source spectrum '^^^^^"'^ for a single GRB is a function of various GRB parameters: 
, T , Td , , ip/ie and e^^. We average over the "measurable" GRB parameters , F , 



and using the procedure described in Ref. j25|] using the same distribution functions for 



these parameters used there |32|. The break energy e^t can be related to total energy in 



photons, £^ (or equivalently to luminosity L^) through the empirical "Amati relation" 33 1 
given by (e^b/lOOkeV) = (3.64±0.04)(,?^/7.9 x 10^^ grg)0-5i±o.oi_ ^^^^ as already mentioned, 

the ratio ^p/ C,e remains as a free parameter. 

What remains to be specified is the GRB rate as a function of redshift, -Rgrb(-2)- Stellar 
core-collapse origin of GRBs as evidenced by CCSN-GRB associations implies that GRB 
rate should follow CCSN rate, Rccsk{z), which is proportional to SFR, Rsf{z). Recently, 
however, an analysis (26] of a reasonably large sample of GRBs with known redshifts from 



the Swift mission 



34| . together with recent accurate determination of the star formation 



history [20|, has given strong indication of a possible enhanced evolution of the GRB rate 
(with redshift) relative to SFR. Taking cue from this we shall allow for a possible effective 
evolutionary factor in the GRB rate relative to SFR and write 



RGRB{z)(x{l + zrRsF{z), 

where a > is a constant, and -Rsf(^) (rate per comoving volume) is taken as 

(1 + for z < 0.97 

Rsf{z) oc { {1 + z) 



(9) 



20 



26| 



(10) 



for 0.97 < z < 4.48 
1 + zy^-^ for 4.48 < z , 

with i?sF(0) = 0.0197 Mq yr^^ Mpc^'^. This SFR including its normalization has been de- 
rived from and is in concordance with recent accurate data on a variety of different indicators 
of SFR in the Universe, and is also in conformity with the experimental upper limit on the 
DSNuB flux giwn by the Supe.Ka,„,oka„de experiment Q. FoUowng the terminology 
introduced in pj| we shall refer to the above SFR as the "concordance model" of SFR. 

The core-collapse supernova rate, Rccsk{z), is related to Rsf{z) through the Initial Mass 
Function (IMF), dn/dm, giving the differential mass distribution of stars at formation. Thus, 



R> 



CCSN 



100 Af0 dn 
8Mq dm 



dm 



Rsf{z) , 



(11) 



where, following standard practice, the IMF is assumed to be epoch (redshift) independent 



l^see, e.g., Il9|). and we have assumed that all stars more massive than ~ 8 Mq undergo 
core-collapse and die on a time scale short compared to Hubble time. Also, our results are 
insensitive to the exact value of the upper cut-off of the IMF (chosen to be at 100 Mq above) 
as long as it is sufficiently large ( >, 30 or so). The SFR (ITUi) assumes an IMF of the 



form 



35(1 . dn/dm oc m "^'^^ for m > 0.5 Mq, and dn/dm (x m ^'^^ for 0.1 Mq ^ m < 0.5 Mq. 



(12) 



With this, the GRB rate can be written in terms of CCSN rate as 

-Rgrb(^) = /grb/ccsn(2;)-Rccsn(^) = /grb/ccsn(0) {1 + z)" -Rccsn(^) , 

where the the normalized core-collapse event rate in the Universe, Rccsn{z), using equations 
(dOD and (dH), is 

^ (1 + z)3-44 for z < 0.97 

12.29 (1 + 2)-°-26 for 0.97 <z< 4.48 (13) 
4.57 X 10^ (1 + z)-'^-^ for 4.48 < z . 



Rccsf^{z) = 2.60 X 10"Vr-^Mpc"^ < 
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The analysis of Ref. 26| seems to indicate the best-fit value of the evolution index a 
appearing in equations ([9]) and f|T2l) to be ~ 1.5, but in this paper we shall keep a as a free 
parameter and study the dependence of our derived upper limits on /grb/ccsn(0) on a. 



III. UPPER LIMITS ON /grb/ccsn 

The DGRBNuB (all flavor) flux calculated from equation ([HD together with equations ^ 
- ([5]), (fT2l) and ( fT3l) with /grb/ccsn (0) = 1 and z^^-^ = 6 (there is negligible contribution 
from z beyond this value), and averaged over the GRB parameters in the manner described 
in the previous section, is shown in Figures [1] and [2] (the superscript "ob" has been dropped 
in these Figures). Figure [1] shows the flux for the equipartition case of ^p/^e = 1 (i-e., 

= = = 1/3 with our choice of C,e = C,b) for five different values of the GRB rate 
evolution index a including the case a = (no evolution), while Figure [2] shows the flux 



for different values of the parameter ^p/^e with a = 1.5, its "best-fit" value from Ref. 26|. 
In both Figures, we also show the current all flavor 90% C. L. upper limit on e^$(ei,) from 



the AMANDA-II experiment [ll| [36|] and also the projected upper limit from the IceCube 
experiment after three years of operation [37], both for an assumed spectrum of the form 
$(e,^) oc e^^. The resulting upper limits on /grb/ccsn(0) obtained by requiring that e^$(ei/) 
not exceed the AMANDA-II limit are shown in Figures [3] and H] as functions of the parameters 
a and ^p/^e, respectively. For comparison, the range of current estimates of the value of the 
ratio /grb/ccsn (0) derived from various astronomical observations 0, Q is also indicated 
in Figures [3] and HI 

It is clear that for a given value of a a higher value of the ratio ^p/^e implies a higher 
predicted level of DGRBNuB flux (see equation ([2])), thus giving more stringent constraint 
on (i.e., a smaller upper-limit value of) /grb/ccsn(0). Similarly, for a given value of ^p/^e, a 
higher value of a implies more GRBs at higher redshifts, again implying a higher predicted 
level of DGRBNuB and consequently more stringent constraint on /grb/ccsn (0). Thus, the 
most conservative limit on /grb/ccsn (0) comes from the case a = and ^p/^e = 1- These 
limits are /grb/ccsn(0) < 5.0 x 10-■^or a = 0, and /grb/ccsn(0) < 1.1 x 10"^ for a = 2. 



For the "best-fit" value of a = 1.5 |26|, we have /grb/ccsn(0) < 1.7 x lO^^ for ^p/^e = 1- 
We also see from Figures [3] and H] that, for a wide range of other values of the parameters 
^p/^e and a, the upper limits on /grb/ccsn(0) derived here from the consideration of high 
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3 4 5 6 7 8 9 
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FIG. 1: The total (all flavor) DGRBNuB flux with /grb/ccsn(0) = 1 for the equipartition case of 
1 for various values of a, the effective evolution index of the cosmic GRB rate relative to 



cosmic star formation rate. The current 
by the AMANDA-II experiment 
after three years of operation 



C. L. upper limit on the diffuse neutrino flux given 



111 ] and the projected upper limit from the IceCube experiment 
37l | are also shown. 



energy diffuse neutrino background are already more restrictive than the current upper limit 
(~ 2.5 X 10^'^) on /grb/ccsn(0) inferred from other astronomical considerations |9l. Il0|. 



At this point it should be mentioned that the AMANDA-II limit we have used above 
actually applies specifically to an assumed diffuse neutrino spectrum of the form $(ej,) oc e^;^. 
The DGRBNuB spectra shown in Figures [1] and [2] are clearly not of this form. Strictly 
speaking, therefore, we should calculate the experimental "AMANDA-II" upper limit for 
our form of the DGRBNuB spectrum and then use that to derive the upper limits on 
/grb/ccsn(0). This can in principle be done by feeding the DGRBNuB spectra calculated 
above to the detailed detector simulation and optimized signal event selection procedures for 
the AMANDA experiment. Clearly, this is beyond our scope in this paper. However, use of 
the AMANDA-II limit in our case here may not be too bad an approximation as a first 
step since, according to the signal event selection criteria of the AMANDA-II experiment 



as explained in Ref. 



it seems reasonable to expect that the dominant contribution to 



the would-be signal events for our spectrum would come from the region around the broad 
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FIG. 2: Same as Fig. (H but for the fixed value of a = 1.5, and various different values of the 
parameter £,p/Ce- 

peak of the spectrum where indeed $ oc e~^, approximately. Thus, while we recognize 
that the upper limits on /grb/ccsn(0) derived here from directly using the AMANDA-II 
limit in our case should be treated with caution, we do not expect significant changes in our 
results (by say more than a factor of few) under a more proper evaluation of the experimental 
"AMANDA-limit" for our spectrum. 

It is interesting to note from Figures [3] and H] that the conservative upper limit on 
/grb/ccsn(0) (obtained with ^p/^e = 1) for the case of a = 1.5, the best-fit value of the 
evolution parameter 



|, is not far above the current estimate of the lower limit on this ratio 
inferred from other considerations. For larger values of ^p/^e the upper limits are even closer 
to the otherwise estimated lower limit on /grb/ccsn(0). This implies that the predicted 
DGRBNuB flux should be detectable by the upcoming detectors such as IceCube which 
will have significantly improved sensitivity over that of AMANDA, unless the estimates of 
/grb/ccsn(0) from direct astronomical observations are gross overestimates (which is pos- 
sible, for example, due to incorrect estimates of the average GRB beaming factor), or that 
the assumption of proton acceleration to ultrahigh energies within GRB jets is invalid, or 
both of these. 

A caveat in the analysis presented above is that it is based on the standard assumption 
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FIG. 3: Upper limits on the fraction /grb/ccsn(0)i obtained by requiring e^$(ej^) for the DGRB- 
NuB to not exceed the AMANDA-II Hmit, are shown as function of the GRB evolution index a 
for various values of the parameter ^-pj^e ■ The shaded region indicates the range of values of 
/grb/ccsn(0) estimated from other astronomical considerations. 

of variability timescales of GRBs on the order of milliseconds, which implies small emis- 
sion regions and consequently large internal targ et photon densities for efficient neutrino 
production through photohadronic processes [38!]. While millisecond timescale variability 
has been seen for many GRBs, this may not always be the case. Efficiency of high energy 
neutrino production in GRBs in the coUapsar model with variability on larger timescales 



has been studied, for example, in Refs. 



24 



391]. Also, neutrino production can be effectively 



quenched in individual GRBs if F, the bulk flow Lorentz factor, is sufficiently large. Clearly, 
more precise determination of the distribution of the bulk flow Lorentz factor and variability 
timescale of the GRBs will be useful in calculating the expected level of the diffuse neutrino 
flux from GRBs more reliably which, together with the results from experiments such as Ice- 
Cube, should be able to place more precise constraints on the fraction of all stellar collapse 
events that give rise to GRBs. 
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FIG. 4: Same as Fig. [3l but as function of the parameter ^p/S^e for different values of the GRB 
evolution index a. 

IV. SUMMARY AND CONCLUSIONS 

In this paper we have attempted to derive upper limits on the fraction /grb/ccsn of all 
stellar core-collapse events that give rise to GRBs, by using the current experimental upper 
limit on the high energy (TeV - PeV) diffuse neutrino background given by the AMANDA- 
II experiment in the South Pole, under the assumption that GRBs are sources of such 
high energy neutrinos. High energy neutrinos are predicted to be produced within GRB 
jets through photopion production by protons and subsequent decay of the charged pions, 
provided protons are accelerated to ultrahigh energies at the internal shocks within GRB jets. 
In our calculation we have allowed for a possible evolution of the cosmic GRB rate relative 
to star formation rate. For a wide range of values of various parameters, the upper limits on 
/grb/ccsn (0) derived here from the AMANDA-II results are already more restrictive than 
the upper limit on this ratio inferred from other astronomical considerations, thus providing 
a useful independent probe of and constraint on the CCSN-GRB connection. The closeness 
of the upper limits on /grb/ccsn (0) derived here (in particular for the case of enhanced 
evolution of the GRB rate relative to the star formation rate at high redshifts) to the 
lower limit on this ratio inferred from various astronomical considerations seems to indicate 
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that the predicted DGRBNuB flux should be detectable by the upcoming detectors such as 
IceCube which will have signiflcantly improved sensitivity over that of AMANDA-II. On the 
other hand, non-detection of the DGRBNuB by the IceCube detector after three years of 
operation, for example, will give more stringent upper limits on /grb/ccsn; but at the same 
time will also imply that either the values of /grb/ccsn inferred from direct astronomical 
observations have been signiflcantly overestimated (which is possible, for example, due to 
incorrect estimates of the average ORB beaming factor) or that the assumption of proton 
acceleration to ultrahigh energies within GRB jets is invalid, or both of these. However, 
more precise determination of the distribution of some of the crucial GRB parameters such 
as the bulk flow Lorentz factor and variability timescale of the GRBs will be needed to 
reliably calculate the expected contribution of the GRBs to the high energy diffuse neutrino 
background, and thereby to determine the upper limits on /grb/ccsn more reliably. To 
conclude, then, the up-coming large volume neutrino telescopes hold immense promise of 
yielding signiflcant information both on the nature of the fundamental physical process of 
particle acceleration in GRB sources as well as on the rate of occurrence of these events in 
the Universe. 

One of us (PB) wishes to thank Nayantara Gupta for helpful clariflcations. 
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